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Abstract 

In calculations of heavy-atom molecules with the shape-consistent Relativistic 
Effective Core Potential (RECP), only valence and some outer-core shells 
are treated explicitly, the shapes of spinors are smoothed in the atomic core 
regions and the small components of four-component spinors are excluded 
from calculations. Therefore, the computational efforts can be dramatically 
reduced. 

However, in the framework of the standard nodeless radially local RECP 
versions, any attempt to extend the space of explicitly treated electrons more 
than some limit does not improve the accuracy of the calculations. The errors 
caused by these (nodeless) RECPs can range up to 2000 cm~^ and more for the 
dissociation and transition energies even for lowest-lying excitations that can 
be unsatisfactory for many applications. Moreover, the direct calculation of 
such properties as electronic densities near heavy nuclei, hyperfine structure, 
and matrix elements of other operators singular on heavy nuclei is impossible 
as a result of the smoothing of the orbitals in the core regions. 

In the present paper, ways to overcome these disadvantages of the RECP 
method are discussed. The developments of the RECP method suggested 
by the authors are studied in many precise calculations of atoms and of the 
TIH, HgH molecules. The technique of nonvariational restoration of elec- 
tronic structure in cores of heavy atoms in molecules is applied to calculation 
of the P, T-odd spin-rotational Hamiltonian parameters including the weak 
interaction terms which break the symmetry over the space inversion (P) and 
time-reversal invariance (T) in the PbF, HgF, BaF, and YbF molecules. 

SHORT NAME: GRECP method: Theory and calculations. 
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(Pseudopotential), Ab initio relativistic method. Electronic structure calcu- 
lation. Molecules with heavy atoms. 

31.15.-hq, 31.20.Di, 71.10.-Fx 

Typeset using REVT^ 



1 



I. INTRODUCTION 



The Dirac- Coulomb (DC) Hamiltonian which is used for calculation of heavy atoms and 
heavy-atom molecules has the form (e = me = h = 1): 

= Eh^(0 + E;^, (1) 

i i>j *J 

where the one-electron Dirac operator is 

= -^c{a■V)+c\(3-l) + Vnuc, 

a, (3 are the 4x4 Dirac matrices, c is the speed of light and Vnuc is the nuclear potential 
which can take account of effect of finite nuclear size etc. 

Higher approximation levels involve modifying the instantaneous Coulomb interaction be- 
tween electrons. These modifications are derived from the Quantum ElectroDynamic (QED) 
theory. The simplest approximation is obtained by including the one-photon exchange. In 
the Coulomb gauge it leads to so-called Dirac-Coulomb-Breit (DCB) Hamiltonian (see, e.g., 

i) 

^DCB ^ ^DC ^ ^ (2) 

where 

Bij{uJij) = }_^[-{ai ■ aj) '—^ + [ai ■ Vi){aj ■ Vj) ] , 

i>j ^ij'^^ij 

and Uij denotes the transition frequency between the electrons i and j. 

A low-frequency expansion of the cosines yields the incomplete Breit interaction Bij{0): 



Bij{0) = -^[di ■ dj + {di ■ fij){dj ■ fij)/rl]/rij 



These terms describe the instantaneous magnetostatic interaction and retardation in the 
electric interaction between electrons. 

In calculations on heavy-atom molecules, the DC and DCB Hamiltonians are usually 
replaced by an effective Hamiltonian 

H^/ = ^[h^^'^^(^J + u^/(^J]+ E (3) 

im iw>jw '^'^wjw 

written only for valence and some outer core electrons denoted by indices and j^', U'^-^ 
is an Relativistic Effective Core Potential (RECP) operator simulating interactions of the 
explicitly treated ( "pseudo- valence" ) electrons with those which are excluded from the RECP 
calculation. In Eq. (|^), h'^'^''^ is the nonrelativistic one-electron Schrodinger operator 

y^Schr _lv^ + V 

'"I — 2 nuc 1 

Contrary to the four- component (relativistic) wave function used in DC(B) calculations, the 
pseudo- wave function in the RECP case can be both two- and one-component. 
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II. FEATURES OF SHAPE-CONSISTENT RECPS 



The RECP approximations with the radially local operator for "shape-consistent" (or 
"norm- conserving" ) pseudoorbitals (pseudospinors) pj-Q are the most widely used in modern 
calculations of molecules with heavy elements. Below we discuss mainly these RECP versions 
and their developments because they provide the largest computational savings for the same 
level of accuracy as compared with other RECP versions. 



A. Advantages 

• Chemically inactive electrons are excluded from RECP calculations. 

• The valence orbitals (spinors) are smoothed in the core regions of heavy atoms to 
generate nodeless pseudoorbitals (pseudospinors). Therefore, the number of the one- 
electron basis functions may be minimized, thus reducing dramatically both the num- 
ber of two-electron integrals and the computational time. 

• The small components of four-component spinors are eliminated and the nonrelativistic 
kinetic energy operator is used. The RECP method allows one to use a well-developed 
nonrelativistic technique of calculation and relativistic effects are taken into account 
with the help of spin-dependent semilocal potentials. As a result, the most part of 
difficulties inherent for the DC(B) molecular calculations can be avoided. 

• In principle, correlations of the explicitly treated electrons with those which are ex- 
cluded from the RECP calculation can be considered within "correlated" RECP ver- 
sions. Reducing the number of explicitly correlated electrons with the help of the 
correlated RECPs is a very promising way to minimize efforts when performing high- 
precision molecular calculations. 



B. Disadvantages 

• By now, different schemes of the RECP generation are suggested which use the radially 
local operator for the effective potential. However, as it was demonstrated in many 
test calculations, they provide the comparable level of accuracy for the same number 
of the explicitly treated electrons. 

It is clear that the explicit inclusion of the outer core electrons into the RECP calcula- 
tion is the way to increase the computational accuracy. However, the extension of the 
space of these electrons more than some limit does not improve the accuracy as is ob- 
tained in all our atomic calculation with RECPs and PPs. Although the errors caused 
by the nodeless RECP approximations can be small enough in competent treatments, 
they still range up to 2000 cm~^ and more for the dissociation and transition energies 
even for lowest-lying excitations that can be unsatisfactory in many applications. 

One can conclude that the scheme of smoothing and the number of the explicitly 
treated electrons are not only responsible for refining the accuracy of the RECP cal- 
culations. 
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• The reliability of the radially local RECP versions is not high for transitions with exci- 
tations from the outer core shells, for the chemistry of transition, rare-earth elements, 
etc. 

• Moreover, the direct calculation of such properties as electronic densities near heavy 
nuclei, hyperfine structure, and matrix elements of other operators singular on heavy 
nuclei is impossible as a result of the smoothing of the orbitals in the core regions of 
heavy elements. 

To overcome these disadvantages of the shape-consistent RECP method, a few develop- 
ments were done by the authors concerning the RECP operator first of all ||^. A detailed 
theoretical analysis of the RECP method is presented in 

C. Historical background of RECP method 

Some of the most important achievements in the development of the effective core po- 
tential (or pseudopotential) theory are listed below: 

• The first ideas and papers devoted to the PseudoPotential (PP) theory of Hellmann 
(1935) [|10[ and Gombas (1935) were published more than 60 years ago. 

• The procedure of treatment of the valence electrons in the frozen sea of core electrons 



was suggested by Fock et al. [12| in 1940. 



The smoothing of orbitals in atomic cores was suggested and the orbital angular mo- 
mentum (Z) dependence of the PP operator for the case of pseudoorbitals smoothed 
in the cores was first emphasized in the pioneering paper of Phillips & Kleinman 
(1959) 



The PP operator with the angular projectors was used first by Heine and Abarenkov 
(1964, 1965) fl^ to generate one-electron model potentials mainly for the solid-state 



calculations. 

Improvement of the PP accuracy when treating several valence electrons explicitly was 



suggested by Weeks & Rice |1T5[ in 1968. 



Goddard III suggested a scheme (1968) [|l^ for the generation of the radially local 



effective potentials on the base of inverting the atomic Hartree-Fock (HF) equations. 

The ab initio radially local effective potentials were efficiently applied by Kahn et al. 
(1976) |]l3 to molecular calculations. 



The nonrelativistic shape-consistent ECP was suggested by Durand & Barthelat (1975) 
[0]. Christiansen et al. (1977) [Q have modified their smoothing scheme. 

The Relativistic ECP method for two-component atomic pseudospinors was proposed 
by Lee et al. (1977) 0. 
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Hafner & Schwarz (1979) [|T8| have split the RECP operator on the spin-independent 
and spin-orbit parts that allowed one to take into account the spin-orbit interaction 
only at the final stage of a calculation thus reducing computational efforts. 



D. Notations 

Below we shall use the following notations for the principal quantum numbers n of 
orbitals: Uf for the inner core (IC) orbitals, Uc for the outer core (OC) orbitals, for 
the valence (V) orbitals, and for the virtual (A) orbitals. The quantum numbers of the 
orbital momentum, total momentum, and its projection will be designated by the /, j, and 
m indices, respectively. We shall write one-electron states as 

for the pseudospinors and as 

for the four-component spinors where 

r = 2j-/, r = (f,a) = (r,fi,a), 

ifnij are the radial parts of the pseudospinors, Pnij and Qnij are the radial parts of the large 
and small components of the Dirac spinors, and Xijm is the two-component spin-angular 
function. 

We shall assign an orbital to the OC or V subspaces using the average radius criterion 
rather than the orbital energy criterion. The difference in such a definition is essential mainly 
for transition metals, lanthanides and actinides in which orbitals with different principal 
quantum numbers may have close orbital energies. For the problems of the RECP generation 
and application, which we are going to discuss, the space criterion is more appropriate. 

III. GRECP OPERATOR 

Following basically the scheme developed by K. Pitzer's group (Lee et al. IQ, Chris- 
tiansen et al. 0), the numerical pseudospinors (pnij{r) are constructed of the large compo- 
nents Pnijir) of the outer core and valence Dirac- Fock (DF) spinors so that the innermost 
pseudospinors of them (for each I and j) are nodeless, the next pseudospinors have one node, 
and so forth: 

~ ( \ - j Pnij{r), r > Rc, 



0,1,. ..,L, j = |/±i 



2 h 
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where L is one more than the highest orbital angular momentum of the IC spinors. The 
leading power 7 in the polynomial is typically chosen to be close to L in order to ensure a 
sufficient ejection of the pseudovalence electrons from the IC region. 

To derive the GRECP components Unij{r), the HF equations are inverted for the V and 
OC pseudospinors (for each / = 0, . . . , L and j = |/ ± ||) so that ipnij are solutions of the 
nonrelativistic-type HF equations in the jj^coupling scheme for a "pseudoatom" with the 
removed IC electrons (Goddard III, 1968) 



nlj\ 



1 (P Id + 

2 dr'^ r dr 2r^ 
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+ £nlj ]^nlj{r) + ^n'nlj^n'lj{r) 



J(r) + K(r) 

r 



(5) 



where Z is the nuclear charge, J and K are the Coulomb and exchange operators calculated 
with the pseudospinors, Enij is the one-electron energy of the corresponding spinor, and en'nij 
are the off-diagonal Lagrange multipliers. 

In general, different potentials Unij are obtained for different pseudospinors ifnij with the 
help of Eq. (|^). It is necessary to construct such an RECP operator which would act on each 
individual pseudospinor as the corresponding potential. This operator may be presented in 
the non-Hermitean (right-handed) form as: 



1+1/2 



u^^ = EE E 



Unlj{r)Pnlj, 



(6) 



n l=0j=\l-l/2\ 



where 



nlj 



\nljm){nljm\, 



i.e., \nljm){nljm\ is the projector on the V'nZjm pseudospinor and the n index runs over the 
OC, V, and A pseudospinor subspaces. The direct use of expression (|]) in calculations is 
impossible because it includes the infinite summations on the n and I indices. 

The Unij potentials contain the contribution from the Coulomb and exchange interac- 
tion with the IC electrons, the contribution from smoothing the spinors, and the contribu- 
tion from relativistic effects as a result of the replacement the Dirac Hamiltonian by the 
nonrelativistic-type Hamiltonian. 

With the help of the property 

{'^nlj\Unlj — Un'ljWn'lj) = 0, (7) 

when neglecting small terms, the operator (^ (Generalized RECP or GRECP operator) can 
be written in the Hermitian form: 
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where 



1=0 j=\l-l/2\ 
L 1+1/2 

Uc l=0j=\l~l/2\ 

L ;+i/2 



EE E 

nc,n'^ 1=0 j=\l-l/2\ 



Unjjir) + Un'Jjir) 



UnMr)]} 



n'Jj, 



(8) 



ml 



is the ordinary spin-angular projector and the first hne in Eq. (||) presents the standard 
radially local RECP operator. It can be shown that the action of this operator on the 
OC, V, and A pseudospinors is equivalent (with the required accuracy) to the action of the 
corresponding potentials on them. 

Note that the nonlocal terms with the projectors on the most important correlating 
functions '^n°'/^j('") localized mainly in the OC and V regions and with the corresponding 
potentials, Un^i^j{r), in principle (i.e. in the cases in which it is essential), can be taken into 
account in the expression for the GRECP operator additionally to the terms with the OC 
projectors Pnjj- 

With the help of the following identities for the P;^ projectors (Hafner & Schwarz, 1979) 



1^ 



P«,i=«±l/2(^, 0-) 



21 + 1 

the GRECP operator can be easily written in the spin-orbit representation. 



(9) 



A. Calculations of Hg and Pb 



The accuracy of the GRECP approximation was tested in correlation structure calcu- 
lations of the Hg and Pb atoms in which we were able to use very flexible basis sets. The 
all-electron Fock-space Relativistic Coupled Cluster (RCC) method has been described in 
and its references. RECP errors in reproducing the all-electron RCC transi- 



review 
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tion energies in Hg for the case of 20 correlated electrons and the equivalent correlation 
[7,9,8,6,7,7] basis sets (for all-electron and RECP calculations, see [^) are presented in 
Table | (compiled from pO|). The highest absolute error in reproducing the transition ener- 
gies with excitation or ionization of single electron is 94 cm~^ for the GRECP 0, 729 cm~^ 
for the RECP of Ross et al. and 1747 cm~^ for the energy-adjusted PP ||2^. The same 
number of electrons, twenty, is explicitly treated in all these RECP versions. The larger 
errors for RECPs |^ and ||22| are mainly due to the neglect of the difference between the 
outer core and valence potentials in these RECP versions (see for more details). 
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Calculations of transition energies in the Pb atom are performed for five lowest-lying 
states in (see Table |l| here) by the four-component and two-component (with 
GRECP) versions of the combined method of the Configuration Interaction (CI) and 
the second-order Many-Body Perturbation Theory (MBPT2) p4|. The 22e-GRECP version 
for Pb was tested and very good agreement of the GRECP/MBPT2/CI calculations with 
the four-component DC/MBPT2/CI results is obtained. 



IV. FROZEN CORE APPROXIMATION FOR OUTER CORE SHELLS 

To perform precise calculations of the chemical and spectroscopic properties, the cor- 
relations should be taken into account not only within the valence regions of heavy atoms 
and heavy-atom molecules but within the core regions and between the valence and core 
electrons as well. In practice, the goal is to achieve a given level of accuracy correlating as 
small a number of electrons as possible, thus reducing the computational effort. However, 
the performance of the RECPs generated for a given number of explicitly treated electrons 
can not always satisfy the accuracy requirements of the correlation structure calculation. 



It can be illustrated on example of Tl p6|. To attain a good accuracy, e.g., within 400 
cm~^ for a group of low-lying states of molecules containing Tl, one should correlate at least 
13 electrons of Tl (i.e., including the 5d shell) and it would be optimal to use the RECPs 
with this number of electrons treated explicitly (i.e., 13e-RECP) like the RECP of Ross 
et al. [21] or our valence RECP version 0. However, all the known nodeless 13e- RECPs 



cannot ensure the above-mentioned accuracy and one should use the RECPs with at least 
21 electrons (e.g., 21e-GRECP for Tl of our group (Q). The 5s and 5p pseudospinors can be 
treated as frozen, providing, nevertheless, the above-declared accuracy. The bp orbitals have 
the orbital energies about four times higher and their average radii are 1.4 times shorter than 
those for the bd orbitals. Moreover, their angular correlations are suppressed as compared 
with those for the bd shell because the most "profitable" polarization states {bd in this 
case) are completely occupied in the lowest-lying states. Therefore, the 5s, 5p states are 
substantially less active in chemical processes. 

One can apply the energy level-shift technique in order to "freeze" the 5s and 5p pseu- 
dospinors P,p6|. Following Huzinaga et al. [^, one should add the matrix elements of 



the SCF field operators (the Coulomb and spin-dependent exchange terms) with these OC 
pseudospinors to the one-electron part of the Hamiltonian together with the level-shift terms 

H Bn,^.ij J2\^cflcf3m){nc^.lcfjm\, (10) 

where Bn^^ij is at least of order |2e„^^;j| of magnitude and Sn^^ij are the orbital energies of 
the OC pseudospinors ipn^jijif) which are frozen. These nonlocal terms are needed in order 
to exclude a collapse of the molecular orbitals to the "frozen" states (the 5si/2, 5pi/2,3/2 
pseudospinors for Tl). 

All terms with the frozen core pseudospinors described here (the Coulomb and exchange 
interactions, the level-shift operator) can easily be presented in the spin-orbit form with 
the help of Eq. (H), as was said above with respect to the GRECP operator (the spin- 
orbit form of the GRECP operator can be found, e.g., in [P|). More importantly, these OC 
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pseudospinors can be frozen (as spinors) in calculations with the spin-orbit basis sets and 
they can be frozen at the stage of calculation of the one-electron matrix elements of the 
Hamiltonian (as is implemented in the MOLGEP code |2^). Thus, any integrals with the 



indices of the frozen spinors are completely excluded already after the integral calculation 
step. 

A. Calculation of TIH 

The Multi-Reference single- and Double-excitation Configuration Interaction (MRD-CI) 



method p9| with the Spin-Orbit (SO) configuration selection scheme was applied in 
the first precise GRECP calculation of molecules. In the GRECP/MRD-CI calculations 
of spectroscopic properties for the 0"'' ground state of TIH (see Table |IT| here), the 
21e/8fs-GRECP for Tl was used (i.e., the 21 electron GRECP |^ with 8 electrons occupying 
the frozen OC pseudospinors 5si/2 and 5^1/2,3/2)- The contracted [4,4,4,3,2] basis set for 



thallium generated in ||26| and the [4,3,1] basis for hydrogen (see ^ttp://www.qche'm.pnpi. 



spb.rv) contracted from the primitive (6,3,1) gaussian basis set of Dunning [30] were used 



Fourteen electrons are correlated in the calculation of spectroscopic constants in TIH 
and, as one can see from Table |ITJ very good agreement is found with the experimental 
data contrary to other known calculations. Close results are obtained very recently in the 
GRECP/RCC calculations of TIH in our joint calculations with group of Prof. Kaldor (Tel- 
Aviv University). The GRECP/RCC and GRECP/MRD-CI calculations of HgH and HgH+ 
are in progress now. 

V. CORRELATIONS WITH INNER CORE ELECTRONS AND BREIT EFFECTS 

On example of the mercury atom PDIjPBI the importance of correlation of the inner core 



electrons (removed from RECP calculations) was investigated with the help of the combined 
MBPT2/CI method [||] and with the RCC method [0. It is obtained that at least 34 
external electrons (occupying the IC 4/ shell, outer core 5s, 5p and 5d shells, valence 6s, 6p 
etc. shells) of Hg should be correlated and the one-electron basis set should contain up to h 
angular momentum functions in order to attain a reliable agreement with the experimental 
data for the low excitation energies within 100-200 cm~^, whereas the errors of the gaussian 
approximation of the GRECP components and the effects of different nuclear models are 
negligible for this accuracy. Otherwise, e.g., for the 20 electron RECPs for Hg known from 
the literature, energies of excitations from the ground state can not be calculated with the 
accuracy better than 200-500 cm~^. 

However, our test calculations show that the main contribution from the correlations 



with the 4/ shell is due to the one-electron correction from the self-energy diagrams ||2| 
and, therefore, this contribution can be taken into account for the 20 electron GRECP at 
the generation stage, i.e. the correlations with the IC 4/ shell can be described by the 
GRECP. The technique of generating the "correlated" GRECP is discussed by the authors 
in [0 . When using such GRECPs in calculations, both required basis sets for the correspond- 
ing atoms and computational efforts can be seriously reduced when performing correlation 
structure calculations. 
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Turning back to the TIH calculation, further improvement of the accuracy can be at- 
tained when correlations with the outer core 4/, 5p and 5s shells of Tl and Breit effects are 
taken into account. We expect that this can be efficiently done in the framework of the 
"correlated" 21e/8fs-GRECP version in which 13 electrons are treated explicitly as in the 
above considered GRECP/MRD-Cl calculation of TIH. 

The first correlated GRECP we have generated for the Hg atom in the framework of 
the MBPT2 method. Despite the considerable improvement in accuracy with this GRECP 
version when calculating excitation energies in Hg, the reliability of the applied generation 
scheme was found to be not very high. A new scheme of generation of correlated GRECPs 
in the framework of the RCC method is in progress now. 

Breit and other QED effects are relatively large for lanthanides and actinides. The 
incorporation of these corrections to calculations with RECPs is discussed in our theoretical 
paper on the GRECP method [^. Some recent test calculations have shown that the Breit 
and other QED corrections for outer core and valence electrons can be efficiently represented 
by a one-electron spin-dependent radially local operator and, thus, can be easily incorporated 
into GRECPs. 

VI. RECPS FOR TRANSITION METALS, LANTHANIDES, AND ACTINIDES 

The disadvantage of the commonly used RECP versions for the transition, rare-earth, 
etc., elements is that they still require the explicit treatment of essentially larger number of 
electrons as compared with other ("normal") elements in order to attain the same compu- 
tational accuracy. 

In paper 0, we added "self-consistent" (SfC) terms to the RECP operator which take 
account of the circumstance that the RECPs generated for the different occupation numbers 
of the outermost d or f shells of these elements are somewhat distinguished. Below we 
will refer to these shells as OuterMost Core (OMC) shells following the partitioning made 
in subsection [11 D[ The significant improvement was attained in reproducing the atomic 
excitation energies in comparison with the conventional RECP calculations. 

Whereas the form of the standard radially local RECP operator accounts for the property 
that the dependence of the corresponding potentials Un^ij{r) on the relaxation of the V 
spinors {fnjj) is weak in the bond making, the SfC correction to this operator takes into 
account the fact of a relatively small relaxation of the OC spinors ipucij caused by the change 
of the OMC shell occupation number in a majority of cases. 

The main features of the SfC RECP constructing for the case of the radially local RECP 

are: 

1. The all-electron DP calculations of two generator states with different occupation 
numbers of the OMC d (/) shell are carried out for an atom under consideration. 
These occupation numbers are designated Ni and N2. 

2. The "generating" RECPs U^^ and U^^ are constructed for these generator states 
employing the standard scheme. The operator for a radially local RECP (incorporating 
the SfC corrections into the GRECP is made similarly) has the form 
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U^ = ^le + f/L» + E E [US{r)-U^Arm,, (11) 



J — I' 2 



where -E^^e is the core energy and Ulj (r) are the radial components of the RECP for 
each pair of quantum numbers / and j. They are derived for the OMC d (/) shell 
occupation number N. 

This operator differs from the standard ("nodeless") RECP operator by the E^^^ 
constant. This value corresponds to the contribution of the core electrons (excluded 
from calculations with the conventional RECP) to the total atomic energy and it is 
usually excluded from the RECP calculations because the £"^^6 o^^y shifts all the 
energy levels by a constant and does not influence on the chemical properties. In the 
case of the SfC RECP, inequality E^^].^ ^ -^^re should be taken into account when 
transitions with the change of the occupation number iV are under consideration. 

3. The one-electron SfC RECP operator is written as 

2 A^i-iVaV 2 ) \ 2 J ' ^ ' 

where B is the adjusting parameter and is the average value (^^ | N | \E') of 
the OMC d (/) shell occupation number operator which in the second quantization 
representation is 



jm 



~ '^l2i^nomalomc.jm ^omclomcjm-, (13) 



where Uomc and lomc are the principal and orbital quantum numbers of the OMC shell. 

At every iteration of the MC SCF, etc., procedures, the value can be determined for 
the wave function ^ without a serious complication of the calculation. 

We should add here that in a high-precision electronic structure calculation of a heavy- 
atom molecule, correlation with the electrons occupying the IC or/and frozen OC states (see 
section |^ can be substantial for calculated properties and dependent from the occupation 



numbers of states of V and other OC electrons as was demonstrated in ||2^ for the OC-OC 



correlation of 5s and hd electrons in Hg. Therefore, in principle, the SfC terms dependent 
on the occupation numbers of the V states can appear in the expression for the SfC RECP 
operator ([T2|) . In this case, the use of the "correlated" RECPs U^^ and U^^ also can be 
required in (|T2p. 



A. SfC GRECP calculations of uranium 

In Table |V| (compiled from |^), one-configuration numerical calculations of transition 
and ionization energies in uranium are presented for our GRECP versions, for the RECP 
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as 



of Ermler et al. and for the energy-adjusted pseudopotential of Kiichle et al. 
compared with the DF energies. 

The GRECP correcting terms described in Section ^ allows one to describe accurately 
the V region of uranium when the 5/ shell occupation number close to that for the generator 
state. The addition of the SfC correction increases the accuracy of reproducing the energies 
of transitions without changing the 5/ shell occupation number in two to eight times and in 
those with changing these occupation numbers in 1.5-3 times. One can see from Table ^ 
that the addition of the quadratic SfC (QSfC) correction (the last term in ([T^)) leads to a 
significant increasing (up to 20 times) of the SfC GRECP accuracy for the transitions with 
the change of the OMC shell occupation number. 



VII. TERM-SPLITTING CORRECTION 



As one can see from Table |V|, the above-discussed SfC corrections give no improvement 
in the term splitting of the ^fli2^d\/2^^\/2 configuration of uranium as compared with the 
generating GRECPs. Analysis of these errors shows that the main contribution (about 90 
%) is due to reproducing the OMC spinors with the help of the nodeless pseudospinors. 

Consider the GRECPs with 24 and 32 explicitly treated electrons for which the OMC 
5/ pseudospinors are nodeless and the Qd ones have radial nodes (for the RECPs with the 
14-electron pseudovalence space, for which these OMC shells are described by the nodeless 
pseudospinors, the errors are greater in general and the analysis is more complicated). 

The simplest way to reduce drastically these errors is to generate and use such GRECP 
versions for which the 5/ shell is characterized by nodal pseudospinors, whereas the 4/ 
pseudospinors are nodeless. To reduce computational efforts, the 4/ shell can be treated 
as frozen in the way described in section However, when approximating the nodal 5/ 
pseudospinors, some additional basis functions and computational efforts will be required. 

If the small magnitude of the OMC shell (5/ here) relaxation is taken into account, there 
is another way out which can be applied for the cases of low-lying states. We suggest to add 
the terms 

■\jTSc _ ^ \xiX2,x-iXi\Xl){x2\ ■ |a;3)(x4| 

X\,X2,XS,X4, 

^ ^ {^XlX2,X^X^ ^XlX^,X^X2^^WXz |Xi)(x2|, (14) 



W Xl,X2,X3 



where Xi = {riomcJomc, jiyf^i) Q) to the GRECP operator. These terms are correcting the 
one- and two-electron integrals with the 5/ pseudospinors (and only these integrals). 

The results of our calculations (Table 0) show that the addition of the Term-Splitting 
(TS) correction (called in by the SO correction) allows one to reduce the most serious 



^ The ji and rrii indices run over all the possible values for a given {nomc, lomc)- In the case of the 
partially occupied OMC 5/ shell, the Xi indices correspond to both OC and virtual pseudospinors 
(i.e., occupied and unoccupied one-electron states in the generator state). 
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errors about 10 times for the term splittings, thus reducing the errors for the transition 
energies between terms to the same order of magnitude as are the errors for transitions 
between the states averaged over the different configurations (Table [TV] ). 



VIII. ELECTRONIC DENSITIES AND SINGULAR OPERATORS IN ATOMIC 

CORES 

To evaluate the matrix elements of the operators being singular in the core regions of 
heavy atoms in a molecule after the RECP calculation of this molecule at the first step, 
the proper shapes of the valence four-component molecular spinors must be restored in the 
atomic core regions at the second step. 

The applicability of the two-step algorithms to the electronic structure calculation of 
heavy-atom molecules is a consequence of the fact that the valence and core electron sub- 
spaces may be represented as the subsystems, interaction between which is described mainly 
by some integrated and not by detailed properties of these subsystems. As a result, the use 
of the GRECP allows one to reproduce accurately the original electronic densities in the 
valence region for excited states as it is demonstrated on U in and on Hg, Tl, Pb, and 
Bi in [^]. Another result is that the electronic structure can be restored in the core regions 
of heavy atoms after the GRECP calculation of a molecule with a good accuracy. 

Both the nonvariational procedures (Pacios & Christiansen, 1985; Titov, 1986) p3|j3^ 



and the variational technique (Titov, 1996) |^5| for the electronic structure restoration after 
the shape-consistent RECP calculation were suggested. 



In [p4| , p6H 39|| , the nonvariational restoration technique was applied to evaluation of the 
hyperfine constants and the parity non-conservation effects in the PbF, HgF, BaF, and YbF 
molecules. The molecular pseudoorbitals in the form of linear combinations of atomic pseu- 
doorbitals were evaluated at the RECP calculation stage. Then, the atomic pseudoorbitals 
were replaced by the unsmoothed four-component DP spinors derived for the same atomic 
configurations as the basis pseudoorbitals. The MO LCAO coefficients were preserved af- 
ter the RECP calculation. Such constructed molecular relativistic spinors were used for 
calculation of the P,T-odd spin-rotational Hamiltonian parameters. 

In a general case, the first-order reduced density matrix, {Dpg}, should be evaluated 
after the molecular RECP calculation on the basis set of atomic pseudospinors whereas the 
matrix elements {^^} of the one-electron operators of hyperfine interaction etc., {A^}, 
should be evaluated in the equivalent basis set of atomic four-component spinors. Then the 
physical values are calculated as 

(A^) = E^.<.- (15) 

pq 

The latest ab initio calculations of the hyperfine, P-odd, and P, T-odd constants 
for the BaF and YbF molecules are performed by our group with the help of the 
GRECP/RASSCF/EO scheme [gS],^ (RASSCF is the Restricted Active Space SCF 
method |^0| and EO is Effective Operator technique based on MBPT2 [p4| , p8| ). This tech- 



nique allows one to take into account polarization and correlation in the valence and core 
regions. The electronic wave function of BaF and YbF for the ground (^Si/2) states is cal- 
culated in the GRECP approximation. The molecular four-component spinors in the core 
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region of heavy atoms are restored in the framework of the non-variational procedure. Core 
polarization/correlation effects are taken into account with the help of the atomic MBPT2 
when constructing the effective valence Hamiltonians for the Ba and Yb atoms. 



A. Calculations of YbF 

The results for the YbF molecule by Mosyagin et al. (1998) are presented in Table [VT| 
compared with other results. For the isotropic hyperfine constant A = (Ay + 2A±)/3, the 
accuracy of our calculation is about 3% as compared to the experimental datum. The dipole 
constant = {A\\ — A±)/3 (which is much smaller in magnitude), though better than in 
all previous calculations known from the literature, is still underestimated by almost 23%. 
Being corrected within a semiempirical approach for a perturbation of 4/-shell in the core of 
Yb due to the bond making, this error is reduced to 8%. Our value for the effective electric 
field on the unpaired electron is 4.9 a.u.=2.5 x lO^'' V cm~^. 



One can see from Table that the values of the Wd constant 

Wdd, = 2f Si/2|i/d|'Si/2), 

where Hd describes interaction of the electron EDM (i.e. Electric Dipole Moment) de with 
the internal molecular field E™°^: 



i/, = 24 n ! ■ E'""', (16) 




from the unrestricted DF calculation of Parpia (1998) ^T|, the most recent semiempirical 
calculation of Kozlov (1997) and our latest GRECP/RASSCF/EO calculation are in a 
very close agreement now. It should be noted that the valence electron contribution to the 
Wd in |^T| is in 7.4% agreement with the corresponding RECP-based calculation of Titov et 
al. (1996) Another recent DF calculation of Quiney et al. (1998) |Q gives the value, 
which is two times smaller. 

It is clear from the calculation that the 4/ electrons in the core of ytterbium should 
be explicitly treated in the planned GRECP calculation of YbF by the relativistic coupled 
cluster method in order to obtain the parameters of P, T-odd spin-rotational Hamiltonian 
with better accuracy. 



B. Variational restoration 

In the variational technique of the restoration [^], the proper behavior of the four- 
component molecular spinors in the core regions of heavy atoms may be restored as an 
expansion on the spherical harmonics inside the sphere with an atomic core radius Rrest- The 
outer parts of spinors are frozen after the RECP calculation of the molecule considered. This 
method enables one to combine the advantages of two well-developed approaches: molecular 
RECP calculations (with gaussian basis sets) and one-center calculations (of atomic type 
with, e.g., numerical functions) in the most optimal way. 
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The most promising application of the two-step method is in a possibihty to "spht" the 
correlation structure calculation of a molecule on two consequent calculations in the valence 
and core regions when the OC and V electrons are explicitly involved into the GRECP 
calculation and the OC and IC space regions are treated at the one-center restoration stage. 

Roughly speaking, the computational efforts for the correlation structure calculations 
in the core and valence regions are added within the two-step technique, whereas in the 
ordinary one-step scheme, they are multiplied by each other. 

The efficiency of the two-step technique was confirmed in the BaF and YbF calculations 
Even these not yet perfect calculations have provided a high accuracy (5-10 %) for 



the P,T-odd spin-rotational Hamiltonian parameters (which can be compared with the 50 
% error for the hyperfine constants in YbF with the correlations in the V region taken into 
account within RASSCF), thus confirming a considerable promise of the two-step scheme 
for the correlation structure calculations. 
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TABLES 



TABLE L Transition energies by the 20e-RCC calculations of Hg with different RECP versions 
in the [7, 9, 8, 6, 7, 7] correlation basis set as compared with all-electron calculations. All values are 
in cm~^. 



htate 


AU-el. 




20 el. 




20 el. 


energy- 


(Leading 


finite 


point 


(jKEL/F 


a 


KEur of 


adjusted 


connguration, 


nucl. 


nucl. 


num. g 


aus. 


Koss et al. 


on ^1 "DTDC 

20 el. rr 


Term) 


Transition energies 






Errors'^ 




















37197 


37260 


-16 


-10 


363 


-666 




39029 


39091 


-1 


4 


378 


-416 




43860 


43925 


22 


28 


437 


619 




55036 


55095 


78 


79 


416 


229 


5d^°6s^^Si/2) ^ 


84159 


84215 


28 


36 


482 


-19 














5di06pi('A/2) 


51656 


51734 


-11 


-9 


428 


-796 


5dl06pl('^3/2) 


60802 


60883 


9 


11 


539 


951 


5d''eSo) 


150636 


150720 


-67 


-58 


729 


112 



^The GRECP from reference with the potentials taken in the numerical and gaussian 
forms, respectively. 
'^The RECP from reference 



'^The PP from reference with the corrected by the factors {2l + l)/2 (M. Dolg, private 
communication) . 

'^In this table, errors were calculated as differences between the transition energies from 
the RECP and all-electron calculations for the same number of correlated electrons and 
equivalent basis sets. The point nuclear model was employed. 
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TABLE II. Absolute errors (AE) of Dirac-Coulomb (DC) 4e-CI and 14e,22e-MBPT2/CI cal- 
culations in reproducing the experimental transition energies (TE) for the lowest-lying states of 
Pb (AE exp.); AE of the 22e-GRECP in reproducing DC calculations (AE DC) in cm~^. Basis 
sets are [3, 5, 3, 2] for CI and [4, 5, 4, 2, 2] for MBPT2. 



4e-Cl 14e-MBPT2/CI 22e-MBPT2/CI 



Conf. 


Term 


J 


Experiment 


DC 


GRECP 


DC 


GRECP 


DC 


GRECP 








TE 


AE exp. 


AE DC 


AE exp. 


AE DC 


AE exp. 


AE DC 




3p 




























3p 


1 


7819 


-807 


(26) 


-535 


(56) 


-393 


(69) 




3p 


2 


10650 


-752 


(36) 


-428 


(62) 


-294 


(77) 




Id 


2 


21458 


-1707 


(46) 


-849 


(89) 


-573 


(120) 




is 





29467 


-1553 


(57) 


-270 


(90) 


-33 


(127) 



TABLE III. GRECP/MRD-CI and other calculations of the spectroscopic constants for the 
ground state of TIH. 





Re 


UJe 




Method 


(A) 




(eV) 


SOCIEX: Tl [8,8,5,2] + H [4,3,1] 








(Rakowitz & Marian, 1997 |4|]) 


1.86 


1386 


2.13 


13e-RECP/S0CI: Tl [4,4,4,1] + H [4,2] 








(DiLabio & Christiansen, 1998 ||) 


1.912 


1341 


1.908 


REP-KRCCSD[T]: Tl [4,5,5,1] + H [3,2] 








(Lee et al, 1998 [H) 


1.910 


1360 


2.02 


21e-REP/KRCCSD(T): Tl [4,5,5,1] + H [3,2] 








(Han et al, 2000 Q) 


1.877 




2.00 


21e/8fs-GRECP/14e-MRD-CI Tl [4,4,4,3,2] + H [4,3,1] 








(Titov et al, 2000 |2|) 


1.870 


1420 


2.049 


Experiment (Grundstrom & Valberg, 1938 [[48(|) 


1.866^ 


1390.7 


2.06 


Experiment (Urban et al, 1989 ^) 


1.8721^ 


1391.3 





^Huber & Herzberg (1979) [|^ have published value 1.87 A which can be obtained from the 

rotational constant Bg. 

'^This value is calculated by us from Bg. 
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TABLE IV. Reproducing the all-electron transition energies (in cm ^) between states of ura- 
nium (averaged over nonrelativistic configurations) in one-configuration calculations with RECPs. 







RECP of 


Enerffv- 




Quadratic 


"Frozen 




DF 


Ermler 


adjusted 


SfC 


SfC 


core 








et al. ||3l| 


PP Q 


GRECP 


GRECP 


if) 


if) 


Num. of el-ns 


All 


14 


32 


24 


24 


24 


24 


Conf. 


TV enerffv 






Absolute error 
























7383 


387 


-498 


-35 


-33 


2 


14 




36159 


332 


130 


4 


6 


3 


16 




13299 


-192 


-154 


-3 


-5 


-1 


-16 




17289 


144 


-621 


-31 


-31 


-1 


-5 




42436 


98 


-188 


-18 


-17 





-5 


5/^6^2 


54892 


-121 


-398 


-14 


-15 


1 


-21 


















5 f^7s2 


16483 


1 76 


788 


-723 





54 


187 


5f^7s2 ^ 


















15132 


-738 


-87 


11 


-11 


-16 


-35 


5f'^7s^7n^ 


15016 


90 


-443 


-37 


-26 


-1 


-2 


5f^7s^ 


38913 


82 


-110 


-37 


-22 


1 


3 


5/^6^2 


34022 


-1287 


-153 


28 


-13 


-26 


-62 


5/^6d^7p^ 


32341 


-794 


-457 


-11 


-23 


-17 


-39 


5/^6(i^ 


53637 


-874 


-245 


-21 


-29 


-17 


-39 


5f^7s^6(i^ ^ 


















3774 


3096 


-748 


-17 


-17 


90 


-96 




















12646 


-441 


-626 


-16 


-15 


-5 





5/^75^6^^ 


42638 


-498 


155 


24 


25 


-5 


1 


5f^7s^6d^ 


10697 


608 


-240 


-10 


-10 


13 


1 


5f^7s^6(i^7n^ 


19319 


390 


-826 


-26 


-26 


6 







45478 


402 


-279 


-13 


-13 


6 





5/^6^3 


54986 


1127 


-581 


-14 


-15 


22 


3 


5/^75^6^^ ^ 
















5f^7s^6(i^ 


29597 


11666 


-1526 


-896 


-104 


466 


48 


5/^7s26(i^ 
















5f^7s'^Q(f7n^ 


18141 


-1367 


-778 


46 


49 


-2 


-2 


5/i7s26(i2 


49158 


-1355 


173 


70 


73 


-3 


-2 


5/i7si6d4 


7584 


1655 


-331 


-39 


-40 


22 


14 


5/^75^6^37^^ 


21154 


779 


-1055 


-11 


-11 


16 


10 


5f7s^6d^ 


48146 


909 


-381 


-13 


-13 


17 


10 


5/^6d^ 


54235 


2810 


-782 


-43 


-45 


42 


27 


5/37s26di ^ 
















5/5 


100840 


430 


1453 


-1860 


22 


105 


291 
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TABLE V. Reproducing the all-electron transition energies between terms for different con- 
figurations of uranium (in cm~^) in one-configuration calculations with GRECPs. 



DF 



Num. of el-ns 

Conf., 

Term 



all 



Trans, 
energy 



SfC 
GRECP 



TS-corr. 

SfC 
GRECP 



GRECP 



24 



24 



32 



Absolute 
error 



TS-corr. 
GRECP 



32 



'5/2'-'"'3/2'''l/2 

J = 1 

J = 2 
J = 3 
J = 4 

J = 5 
J = 6 

^/5/2^/7/2'^^l/2 

J = 2 
J = 3 
J = 4 
J = 5 
J = 6 
J = 7 
J = 8 



18643 
9729 
7813 
6759 

-7918 
-10695 



J= 1 



4392 
2843 
3477 
2805 
4631 
-5951 
-5450 



136 
-117 
-140 

-547 
-509 



175 
113 
133 
105 
169 
-275 
-258 



-97 
-56 
-44 
-36 
42 
76 



-52 
-33 

-40 
-32 
-56 
48 
42 



142 
146 
-69 
-94 

-485 
-464 



150 
97 
113 
88 
138 
-253 
-239 



-14 
-13 
-8 
-5 
-11 
-31 



-39 
-24 
-30 
-23 
-42 
28 
24 



5/5/2 6^3/2 '''•^1/2 
J= 1 

J = 2 

J = 3 
J = 4 
J = 5 
J = 6 



J 



: 

-19078 
-15304 
-23607 
-25984 
-32497 
-39551 



-487 
-331 

-682 
-715 

-848 
-785 



-7 
6 
-8 
24 
19 
82 



-461 
-326 
-644 
-689 
-818 
-788 



-64 
-45 
-85 
-74 
-94 
-62 
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TABLE VI. Parameters of the spin-rotational Hamiltonian for YbF. 





A 






Wa 


Ws 


Method 


(MHz) 


(MHz) 


(1025 Ji^) 

^ ecrri'' 


(Hz) 


(kHz) 


Semiemp. (Kozlov et al, 1994 ^) 






-1.5 


730 


-48 


Semiemp. (Kozlov, 1997 Q) 






-1.26 




-43 


GRECP/SCF (Titov et al, 1996 |3) 


4932 


59 


-0.91 


484 


-33 


GRECP/RASSCF 


4854 


60 


-0.91 


486 


-33 


DHF (Quiney et al, 1998 [||) 


5918 


35 


-0.31 


163 


-11 


DHF+CP 


7865 


60 


-0.60 


310 


-21 


DHF (rescaled) 






-0.62 


326 


-22 


Unrestricted DF (Parpia, 1998 Q) 






-1.203 




-22 


DF (unpaired elect- 












ron contribution) 






— u.yoz 






GRECP/RASSCF/EO 












(Mosyagin et al, 1998 H) 


7842 


79 


-1.206 


634 




GRECP/RASSCF/EO 












(with 4/-hole correction) 


7839 


94 


-1.206 


634 




Experiment (Knight et al, 1970 [||] 


7617 


102 
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